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Statistical  Treatment  of  Crack  Propagation  Data 


R.  A.  Heller1 


Summary 


Crack  propagation  data  are  available  for  two  different  particulate  composite  mate¬ 
rials.  A  larger  set  of  38  observations  is  used  to  draw  statistical  inference  for  the  second 
set  consisting  of  5  data.  The  Paris  crack  propagation  relation,  whose  two  parameters 
are  functionally  related,  is  used. 

Introduction 


Two  sets  of  data,  on  particulate  composite  materials  containing  hard  particles  in 
a  rubbery  matrb^Jiave  been  examined;  one  consisting  of  38  observations  and  a  smaller 
one  of  5  observations.  It  is  desired  to  obtain  statistical  inference  for  the  smaller  data 


set,  including  confidence  limits,  with  the  aid  of  the  larger  group  of  observations^Tt 
has  been  postulated  that  the  rate  of  crack-propagation  follows  the  Par-is"felationship 
and  that  the  two  parameters  of  thisrule  are  functionally  related^  It  has  also  been 
indicated  that  the  parameters  are ^log)  normally  and  nonnplly distributed.  The  two 
material^iesignated  as  MM  for  the  larger  set  and  MX ^oftie  smaller,  though  differenty^- 
have  similar  characteristics.  Furthermore^!  M  was  tested  at  a  cross-head  speed  of  0.1 
in/min  and  MX  at  0.2  in/min.  ^ 


According  to  the  Paris  rule,  the  rate  of  crack  propagation 

^  =  Mf2  =  d  .  (1) 

dt 

where  Ci  and  C2  are  correlated  parameters  and  Ki  is  the  stress  intensity  factor  in 
psi  Via- 

Taking  logarithms  of  both  sides 

log  a  =  log  Ci  -r  C2  log  Ki  (2) 

log  Ci  as  well  as  C2  are  assumed  to  be  normally  distributed  random  variables  (C 1  is 
consequently  log  normally  distributed). 

As  a  result,  log  a,  the  sum  of  two  normally  distributed  variables,  is  also  normal 
with  a  mean  of 

log  d  =  log  Ci  4-  C2  log  Ki  (3) 

and  standard  deviation 

^ogd^kfogCx+log  2*V<4]1/2  (4) 

The  relationship  between  log  Ci  and  C2  is  linear 

log  Ci  =  A  —  BC2  (5) 
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or  in  terms  of  standardized  variables 


(C2  -  C2)/crc2  =  -(log  Cl  -  log  Ci)/<7)og  Cj  (6) 

that  is 

log  A  =-^£L(c2-C2)  +  hfC^  (7) 

<7C2 

Substituting  Eq.  7  into  Eq.  2 

log  a  =  c2  ( log  Kj  =  (8) 

V  &C2  J  &c2 

is  obtained.  Therefore,  the  variance 

crfog  d  =  Wc2  log  Ki  -  criog  cj2  (9) 


Data  Analysis 


•JLO.c  £, 


The  two  sets  of  data:  38  observations  of  log  Ci  and  C2  for  MM,  at  a  cross-head 
speed  of  0.1  in/min  and  5  observations  for  MX,  at  a  cross-head  speed  of  0.2  in/min  are 
listed  in  Table  (^together  with  their  sample  means  and  sample  standard  deviations. 
The  crack  propa^a&on-rates,  log  a  for  Ki  =  40  psi  Vm  are  also  listed  in  the  table. 


TableQ)  C2,  log  C\  and  log  a 


O  OajL  ch 

■»  f 

-  r\  UrA-'.r  . 

44vcy*  (COA'  o,  A 


/Vr^b'  c.  » 


Test  No. 

c2 

log  Ci 

log  a 

103 

2.42161 

-8.91898 

-5.03956 

105 

1.76926 

-6.7678 

-3.93344 

106 

1.507719 

-5.85985 

-3.4449 

107 

1.56330 6 

-6.09666 

-3.59224 

108 

1.996521 

-7.48854 

-4.29011 

109 

2.326644 

-8.55145 

-4.82417 

110 

1.718226 

-6.53345 

-3.78085 

111 

1.882954 

-7.08362 

-4.06712 

112 

2.212894 

-8.21514 

-4.67008 

113 

1.885258 

-7.15843 

-4.13824 

114 

2.286922 

-8.46155 

-4.7979 

115 

1.961836 

-7.40744 

-5.80544 

117 

1.758835 

-6.6976 

-3.87995 

118 

2.358755 

-8.72234 

-4.94362 

119 

1.403963 

-5.53812 

-3.28897 

120 

2.205351 

-8.19052 

-4.65755 

121 

1.981772 

-7.43898 

-4.26416 

122 

1.645804 

-6.34225 

-3.70568 

123 

2.259029 

-8.3731 

-4.75414 

124 

1.22662 

-4.96987 

-2.78909 

•  125 

1.867377 

-7.09606 

-4.10452 

128 

1.787407 

-6.82663 

-3.96321 

129 

1.156759 

-4.6236 

-2.77048 

130 

1.961972 

-7.344 

-4.20092 

131 

1.551135 

-6.02706 

-3.54124 

Test  No. 

C2  * 

log  Ci 

log  a 

132 

1.527031 

-5.86331 

-3.41701 

133 

1.827629 

-6.91117 

-3.98331 

134 

1.608571 

-6.17429 

-3.59736 

136 

2.094406 

-7.88458 

-4.52935 

137 

1.169022 

-4.7417 

-2.87893 

138 

1.381992 

-5.43623 

-3.22227 

139 

1.273468 

-5.06151 

-3.02142 

140 

2.019887 

-7.53751 

-4.30165 

141 

1.944194 

-7.29269 

-4.17809 

142 

2.516193 

-9.19896 

-5.16802 

143 

1.994757 

-7.45045 

-4.25485 

144 

1.735048 

-6.99552 

-4.21597 

145 

2.349725 

-8.59174 

-4.82748 

Mean 

1.845786 

-6.99665 

-4.07457 

Std 

'  0.365129 

1.214061 

0.70221 

G13L 

2.857 

-7.554 

-2.97709 

G13R 

4.255 

-10.287 

-3.47049 

G23 

2.775 

-7.359 

-2.91345 

G19L 

3.881 

-9.581 

-3.57748 

G19R 

4.461 

-10.724 

-3.36364 

Mean 

3.6458 

-9.101 

-3.26043 

Std 

0.786049 

1.557153 

0.29832 

ptr'i  Ir'^y  ^ 


■ - The  dat^72  vsriog  Ci,  axe  plotted  for  the  two  materials,  MM  and  MX,  in 

Fig.  1  and  indicate  the  validity  of  linear  relationships  with  correlation  coefficients  of 
p  =  -0.99836  and  -0.9999^respectively. 


Fig.  1.  Log  Ci  vs.  Gi  for  a)  MM  and  b)  MX 

It  is  seen  that  the  range  of  values  are  quite  different.  For  MM,  C2  varies  between 
1.16  and  2.52  with  corresponding  log  C\  values  of  -4.6  and  -9.2.  For  MX,  C2  varies 
from  2.8  to  4.46  with  log  Ci  from  -7.55  to  -10.72. 


The  differences  axe  attributable  to  the  two  different  materials  rather  than  to  the 
changed  cross-head  speeds.  Though  the  cross-head  speed  fpr  MX  is  double  that  for 
MM,  the  mode  of  cracking  is  similar.  It  is  expected  that  only  at  much  greater  cross- 
heads^speeds  would  cracking  change  to  a  brittle  fracture. 

Utilizing  Eqs.  3  and  4,  log  a  has  been  calculated  for  each  datum  for  a  stress- 
intensity  factor  of  Kj  =  40  psi\/in  (Table  1). 

The  data  were  plotted  on  normal  probability  paper  in  Fig.  2.  Though  all  datum 
points  fit  the  normal  distribution  reasonably  well,  all  MX  values  axe  higher  than 
log  a  for  MM. 


In  order  to  eliminate  the  influence  of  the  material  differences,  the  log  a  values  were 
standardized  using  the  appropriate  means  and  standard  deviations  for  each  material 
(Table  1).  The  data  were  then  arranged  in  increasing  order  and  were  plotted  on  normal 
probability  paper  (Fig.  3). 


It  is  apparent  that,  in  this  standardized  fornj^oth  sets  are  essentially  normally 
distributed  and  that  they  belong  to  the  same  populatiohQthe  MX  data  are  interspersed 
with  MM  values.  ( 

4k;.  c  c.->i4k_ 


Ac 


j. 


'covoo  or  o l  «•'  ) 

(dw\) 


Confidence  Limits 


Confidence  limits  may  be  established  for  the  normally  distributed  data.  The 
meaning  of  such  limits  is  that  the  true  population  mean  will  lie  between  these  limits 
with  a  given  probability  (1  —  a).  The  width  of  the  confidence  interval  is  a  function 
of  the  number  of  observations,  n.  The  greater  the  value  of  n,  the  narrower^j^the 
confidence  band. 


For  small  numbers  of  observations,  confidence  intervals  axe  based  on(§tudentJi  t 
distribution  [1],  For  a  variabl^X,  with  sample  mean,  ^and  standard  deviation,  S, 
the  population  mea^/xi_a/^vill  lie  between  the  following  confidence  limits 

[(x  ta/2)  n  —  1  S/y/n)  <  fj*a- 1  <  (JY  +  Ta/2j  71  ~  1  S/y/nj]  (10) 


where  n  is  the  sample  size  and  fa/2,n-i,  n—  1  is  the  tabulated  value  of  the  t  distribution 
with  a/2  and  number  of  degrees  of  freedom  /  =  n  —  1. 


In  the  case  of  log  a,  for  MM ,  with  a  confidence  1  -  a  of  .99  log  a  =  4.07457, 
Siog  d  =  .702212  and  n  =  38,  a  =  .01,  a/2  =  .005,  /  =  37,  p  =  1  -  a/2  =  .995. 

t. 995,37  =  2.72  (by  interpolation) 


-4.07457  -  2.72 


/.70221V 
V  v/38  ). 


^  og  d(l  — q)  < 


—4.07457  4-  2.72 


/. 702212 
V  \/38 


(11) 


Therefore 


-4.38442  <  /ilos  d(1_a)  <  -3.76472 


(12) 


a  difference  of  0.62. 

The  true  mean  lies  between  -4.38442  and  -3.76472  with  a  confidence  of  0.99. 

— -  In  contrast^for  MX  with  log  a  =  —3.26042,  Siog  d  —  .298318,  n  =  5,  for  a 
confidence  of  1  —  a  =  .99,  /  =  n  —  1  =  4,  ta/2,n-i  =  t. 995,4  =  4.604 

(-3.26043  - 4.604  ~29^18^  <  Mlog  d(1_a)  <  (-3.26043  +  4.604  8  ) 

(13) 


(-3.87466  <  /rlog  d(1_o)  <  -2.64620)  (14) 

with  a  difference  of  1.228. 

The  confidence  interval  for  MX  is  much  greater  than  for  MM,  because  the  former 
set  has  a  small  sample  size. 


To  improve  on 


normal  distribution, J7 ig^4/ 
of  unity,  Eq.  11  is  modified 


Lfidence  interval  for  the  smaller  sample,  the  standardized 
^nay  be  used.  With  a  mean  of  zero  and  standard  deviation 
as 


[(0-fa/2,n-il/\/«)  </iQ_i  <  (0+ta/2,n_il/Vn)] 


(15) 


Again  for  1  —  a  =  .99,  with  n  =  43,  t. 995,42  =  2.6948. 


For  the  standardized  variables 
[(0  -  2.6948  1/Vuj  <  n. 99  <  (0  +  2.6948  1/V44j] 
or 


-0.4063  <  /z.99  <  .4063 


(16) 


X/  <  /x.99  <  Xu 

Performing  the  standardization  backwards 
(xi  x) f Sx  —  (xu  X)/Sx  = 


Xwt  "/UL-. 

*'  ^  s  V\Du_- 


cLev  ■  a.'V".  oti 


Where  £  and  w  indicate  the  lower  and  upper  confidence  levels,  x  and  Sx  are  the  mean 

and  standard  deviat  ion  of  log  a  for  MX 
\~S 

Xt  =  log  ai  =  -.29832  x  .4063  -  3.26043  =  -3.3816 


Xu  =  log  du  =  .29832  x  .4063  -  3.26043  =  -3.1392 


Therefore  the  mean  log  a  for  MX  with  a  confidence  of  .99  lies  between 


3.3816  <  /x.99  <  -3.1392 


with  an  interval  width  of  .242  instead  of  1.228  based  on  the  original,  data. 

The  validity  of  the  suggested  calculation  hinges  on  the  assumption  of  normality 
for  both  sets  of  data  and  that  crack  propagation  is  not  significantly  different  in  the 
two  materials. 
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